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ABSTRACT 


iieimaeepaeccrns with the appearance of cyclonic motion 
and sharp thermal fronts frequently are seen in satellite 
teetmases otf the California coast near Pt. Sur. These 
thermal patterns are associated with distinctly structured 
upwelling systems. In nine seasonal cruises since December 
1978, nutrient fronts are strongly correlated with thermal 
fronts. However, an upwelling which is well defined by the 
Saeecllite detected thermal front, may not be nutrient rich. 
A strong inverse linear correlation between nutrients and 
temperature though does exist for upwelling within an early 
stage of development. It is feasible, by utilizing limited 
Mims eeeddta, to inter nutrrent distribution within satellite 


detected, surface thermal patterns associated with upwelling. 
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I. INTRODUCTION 


tiem ihesis 0S a study of the linear correlation of sat- 
ellite detected sea surface thermal patterns to nitrate and 
phosphate levels in a coastal upwelling region. In part, it 
considers whether the chemical mesoscale in upwelling systems 
may be inferred from satellite IR detection of surface thermal 
patterns. 

The Advanced Very High Resolution Radiometer (AVHRR) on 
the TIROS (Television Infrared Observation Satellites) sat- 
ellite series can be used to locate sea surface temperature 
patterns which identify upwelling regions for in situ auto- 
mated biochemical and thermal sampling and analyses. Surface 
ampenatune and nutrient maps can be generated using in situ 
measurements and inference re temperature based on synoptic 
satellite IR imagery and nutrient to temperature correlations. 

Previous studies which used the same methods to detect and 
describe the thermal features (e.g., Traganza, Nestor, and 
McDonald, 1980) showed that some of these features resulted 
from upwelling pulses and produced nutrient structures of the 
same scale. These highly structured upwelling systems had 
Sharp thermal and chemical gradients or fronts. The recurrent 
formation of a "cyclonic upwelling system" off Pt. Sur, Cali- 
fornia [Traganza, Conrad, and Breaker, 1980], offers a unique 


opportunity to investigate the significance of frontal systems in 
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regions of upwelling. Significantly higher nutrient concen- 
trations observed in upwelling systems can enhance biological 
Setivity in the proximity of the fronts. 

The ability to infer the nutrient concentrations associ- 
ated with an upwelling system, identified in satellite IR 
imagery, might help in explaining the relationships between 
the chemical mesoscale and pelagic ecosystem in the region. 
The occurrence of microplanktonic blooms adjacent to sharp 
thermal and chemical gradients has been documented in many 
coastal upwelling regions around the world [Barber, 1980; 
Blasco, 1980; Boyd and Smith, 1980; Traganza, Conrad, and 
Breaker, 1980]. Because phytoplankton are the primary pro- 
ducers in the pelagic ecosystem, their productivity supports 
the foodweb of which man harvests primarily the higher level 
heterotrophs. The magnitudes of nutrient concentrations and 
nutrient fluxes in upwelling systems could provide clues to 
the biological productivity in such areas. The associated 
thermo-chemical gradients or fronts could become the foci 
for fishing efforts. With the potential for biological pro- 
duction known and the growth rate predicted for an ecosystem, 
fishing fleets could be effectively dispatched to optimize 
yields. 

-While thermal patterns may correlate well with satellite 
observations, correlations for chemical patterns with IR 
imagery are by inference and may be less simple, tending to 


vary as a function of the initial conditions, boundary 
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conditions, biological processes, and the dynamic processes 
acting over time. An understanding of all these processes 
may be necessary to translate satellite thermal patterns 
into nutrient levels. 

The objectives of this thesis are: 

Molomtmvestigate the linear relationship between in 
situ values of nutrients (nitrate and phosphate) and temper- 
ature associated with coastal upwelling occurring under 
varying nutrient and temperature characteristics of the 
ambient water mass and stages of development. 

2. To map observed values of temperature, nitrate, 
phosphate, and nutrient ratio (nitrate/phosphate). 

3. To compare the magnitudes and gradients of tempera- 
ture and nutrients associated with the oceanic front in 
different types of upwelling in different stages of develop- 
ment. 

4. To examine the relationship between the surface 
temperatures and nutrient concentrations in recently upwelled 
water and to estimate the depth from which it advected. 

5. To test the accuracy of a satellite-aided prediction 
in which the source water temperature and nutrient character- 
istics are used in a regression equation to hindcast surface 


nutrient concentrations. 
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II. METHODS 


A. STUDY AREA 

Micmsoulcdyedtcan(rag. 1) extends ca. 30 km north of Pt. 
pureand ca. 60 km to the south. It is bounded by the 122° 
30' W meridian and the Central California coastline stretch- 


ing from Monterey to Cape San Martin. 


bee CRUISES 

Three cruises were conducted aboard the Naval Postgraduate 
School's research vessel R/V ACANIA, 27 and 28 September 1979; 
29 and 30 November 1970; and 10 and 11 June 1980. Positions 
were determined by LORAN-C which has good coverage in this 
region. Running plots of temperature and nitrate and phosphate 
concentrations were maintained and used to ensure coverage of 
features of interest. Particular attention was given to the 
location of thermal gradients observed in satellite images 


received prior to each cruise. 


C. SURFACE SEAWATER SAMPLING SYSTEM 

Seawater was continuously sampled from a depth of 2.5m. 
The September and November cruises pumped seawater from a 
keel intake via the ship's pump to the shipboard laboratory. 
The June cruise used an overboard nylon hose intake at 2.5m, 
which was attached to a deck-mounted Wilden, M-4, air oper- 


ated diaphragm pump (Wilden Pump and Engineering Company, 
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Fig. 1. Bathymetry of the study area with grid superimposed. 
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Colton, California 92324). An air debubbler was used in 
line between the pump and measuring equipment [Technicon 
Autoanalyzer, AA-II and Turner III Fluorometer; see Bronsink, 


1980]. 


D. VERTICAL PROFILE SAMPLING 

Seven Nansen casts were conducted also on the November 
cruise. Seawater samples were collected at 25, 50, 100, 200, 
300, and 400 meters for analysis of temperature, salinity, 


and nitrate. 


fee LEMPERATURE 

A thermistor sensed the temperature of the seawater 
entering the keel intake. A strip chart recorder continu- 
ously plotted it. The equipment was calibrated and monitored 
every 30 minutes by bucket thermometer readings of the sea 
surface. The mean difference, T(bucket) - T(thermistor), of 
0.1°C (std. dev. +0.2°C) was determined for the temperature. 
Sippican expendable bathythermograph (XBT) probes were re- 
Peascamcvem, 20 Minutes, or ca. every 5 km, during the Sep- 
tember cruise. The September 1979 and June 1980 cruises 
dropped XBT probes only in areas of major interest. A spac- 
ing of ca. 5 km was maintained along the major axis of the 
feature and during front transits. For calibration, sea 
surface temperature was measured by bucket thermometer 


Simultaneously with the launch of each XBT probe. The mean 
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difference, T(bucket) - T(probe), between surface measure- 


ments was 0.0°C (std. dev. +0.2°C). 


F. NUTRIENTS 

Surface (2.5m) nutrient concentrations were analyzed 
every two minutes according to the Technicon Industrial Method 
175-72-WM, 177-72-WM [Anonymous, 1973] and 100-70-WM [Anon- 
Peonlo7o)|. Nitrate here may include traces of in situ 
nitrite, since the nitrate is reduced to nitrite before meas- 
urement. However, according to Paulson [1972] there is little 


or no nitrite at the surface sampling depth in this area. 


G. LINEAR REGRESSION AND CORRELATION COMPUTATIONS 

Regression lines were generated for nitrate versus 
phosphate, nitrate versus temperature, phosphate versus tem- 
perature, and nutrient ratio (nitrate/phosphate) versus tem- 
perature. The linear regression analysis used the NORLSQ 
library subroutine from the W.R. Church Computer Center at 
the Naval Postgraduate School. 

The population correlation coefficient was derived by 


the equation: 


NLxX5y; z (Zx5) (Zy;) 


we See ee 
({ntx? - (£x,)?](nty? - (Zy,)?])? 


When a "none detected" condition occurred due to concentra- 
tions below the sensitivity of the instrument, the nutrient 


and temperature data were not used in the calculations. 
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Heme ORTELLITE IMAGERY 

The AVHRR has a temperature resolution of .5°C within 
the IR (10.5 um to 12.5 um) spectral window and a spatial 
resolution of 1.1 km. The thermal energy measured by the 
AVHRR is primarily a function of radiation emitted from the 
upper millimeters of the sea surface, clouds, and atmospheric 
gases. The greatest errors in estimating sea surface temper- 
ature are induced by atmospheric attenuation of the thermal 
energy due to absorption and scattering. The effect of 
absorption by moisture in the atmosphere presents the most 
severe problem because of the magnitude of absorption and the 
variability of the moisture field. At mid-latitudes over the 
ocean, this attenuation may result in measurements that are 
2 to 6°C lower than actual sea surface skin temperatures 
[Maul and Sidran, 1973]. 

The AVHRR has a temperature response from ca. -90 to 60°C. 
The measured thermal energy is normally displayed as gray 
tones on photographic film to produce images interpretable 
as Clouds, land, and water. The colder temperatures within 
these images are assigned lighter shades of gray. Image 
enhancement reduced the range of temperatures to those rep- 
resentative of the sea surface, while retaining all the gray 
tones [see Methods: Satellite Assisted Sampling Strategy]. 
In the enhancement of the TIROS-N images (Plates 1, 2, and 
Pmcmceiperature range representative of the in situ range 


of values was assigned to shades of gray with white equal 
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momcaee Ceand black Ga. 22°C. Studying satellite images 
alone is subjective because of the inability to assign tem- 
perature values representative of the sea surface to each 
gray shade. However, the relative temperature values and 
Mecemal gradients are more representative of in situ measure- 
ments within the same image than are absolute values. Because 
of the atmospheric effects on the thermal energy reaching the 
satellite, these measurements can best be described as 


Saoparent temperatures." 


eon LELLITE ASSISTED SAMPLING STRATEGY 

Approximately one week prior to each scheduled departure, 
Mr. Breaker, the staff oceanographer at the National Environ- 
mental Satellite Service (NESS), in Redwood City, California, 
would be alerted. He closely monitored images of TIROS-N 
satellite series for oceanic features of interest in the 
study area. The images were computer enhanced by him to 
better delineate sea surface thermal patterns. All images 
had the non-linear distortion associated with the earth's 
curvature and rotation removed. The day prior to the cruise 
and sometimes during the cruise, Mr. Breaker sent updated 
satellite information regarding the approximate location, 
center, size, and orientation of the feature. When possible, 
the most useful and most recent images were delivered prior 
to sailing. All proved invaluable in locating features and 


planning sampling strategy. 
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J. SATELLITE ASSISTED MAPPING OF SURFACE TEMPERATURE 

Two objective temperature fields were drawn by the CONISD 
Gompucer program; one field was formed from only in situ 
temperature, the other included data from satellite imagery 
ase well. 

The size of the study area was small enough that the 
spherical distortion attributed to the earth's curvature was 
assumed to be negligible. Therefore, the author constructed 
a grid on the assumed planar surface (see Appendix C) with 
the head of the Sur submarine canyon as its origin (Fig. 1). 
Each in situ data point had cartesian coordinates assigned. 
Using this annotated data, the CONISD library subroutine from 
the W.R. Church Computer Center drew surface contours. Be- 
cause the contouring computations involved a triangulation 
process between data points nearest the limits of the grid 
and the study area boundaries, it was necessary to assign 
ocean and near coast boundary values that were realistic to 
generate a surface map which well represented the external 
limits of the surface feature. These values were inferred 
from the thermal patterns visible in satellite IR images. 

Enhanced satellite images taken during or nearest the 
period of in situ observation and having a cloud-free view 
of the study area were used. These images were taken to NESS 
to geographically locate the surface thermal patterns associ- 
ated with the upwelling system on a transverse mercator 


projection of the study area. The author outlined shifts in 
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Memory cones, Visible to the naked eye. These lines rep- 
resented a subjectively determined change in the gray scale. 
A uniform distribution of atmospheric moisture over the small 
area of the surface feature was assumed. The change in gray 
scale could therefore be related to synoptic variations in 
apparent surface temperature. Mr. Breaker, with his exper- 
tise in using the Baush and Lomb Zoom Transfer Scope, trans- 
cribed the apparent thermal patterns onto the transverse 
mercator projection. These patterns were then transferred 
by the author onto the grid developed for the study area. 
Using a surface temperature map generated only from in 
Smaiemeasurements, a plot of the ship’s track (to know where 
and when in situ measurements were taken) and the synoptic 
thermal patterns inferred from the satellite image of the 
feature, the thermal patterns within the area of concentrated 
Buesttu Measurements were assigned temperature values equiva- 
lent to those on the surface temperature map. The satellite 
detected and in situ thermal patterns matched well, within 
+.5°C. Where these thermal patterns approached the coast, 
those values inferred above were applied to the near coast 
boundary. For objective (machine) contouring the oceanic 
temperature field beyond 6 km from the seaward outline of the 
surface thermal feature was arbitrarily assigned a uniform 
value. The oceanic boundary formed a ee at a constant 
distance (ca. 6 km) seaward of the surface thermal feature. 


The warmest in situ measurement made in this oceanic region 
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was assigned to the oceanic boundary. With the inferred 
near coast and ocean boundary values and in situ data, 
another surface temperature map was generated by the CONISD 
subroutine. Hand smoothing was applied to reduce the dis- 
tortion thought to be caused by advection or other dynamic 
processes over the time period of the study. These temper- 


ature surface maps appear in Figures 10, 20, and 34. 


K. SURFACE NUTRIENT MAPS 

The sea surface nutrient maps (Figs. 11, 12, 13, 21, 34, 
35, and 36) were generated in similar fashion to the sea sur- 
face temperature map. Each in situ data point was assigned 
cartesian coordinates on the study area grid. The near coast 
and ocean boundary values for temperature were put into the 
regression equations for nitrate/temperature, phosphate/ 
temperature, and nutrient ratio/temperature. The respective 
nitrate, phosphate, and nutrient ratio output were assigned 
to the coordinates of the input temperature. Using this 
omotatccuin situ and temperature inferred data, the CONISD 
library subroutine drew surface contours. Hand smoothing 
was applied to reduce the distortion thought to be caused by 
advection or other dynamic processes over the time period of 


the study. 


L. DEFINITION OF FRONT AND CALCULATION OF ITS GRADIENTS 
For the purpose of this thesis, the oceanic front is 


defined as the area where the thermal and chemical gradients 
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reach their maximum change per kilometer [see Discussion: 
itemnelatitonsnip Between Source Water Characteristics and 
the Oceanic Front]. In order to locate the oceanic front 
and assess its gradients, the following procedures were 
Samried Out. 

The incremental change per kilometer of nutrient concen- 
trations and temperature was determined from the in situ 
values observed along the ship's track. The elapsed distance 
midway between the magnitudes used to compute this spatial 
gradient was assigned as the position on the ship's track at 
mmeehmune Ghange occurred (Figs. 14, 22, and 38). This tech- 
nique was used to take advantage of the high spatial resolu- 
tion (better than 0.6 km) along the ship's track which showed 
microstructure fronts within the feature as well as the 
Seeanic front. 

The magnitude of the gradient of the front was computed 
after correcting for the angle at which the ship transited 
the front. It was assumed that the front parallels isolines 
of temperature and nutriend concentrations on the machine- 
contoured charts. The angle of incidence the ship's track 
intersected these isolines at the vicinity of the front was 
measured using surface nutrient and temperature maps. Apply- 
ing the law of sines, the distance along the normal was found 


and used to compute the magnitude of the gradient. 
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Dit RESULTS 


PeEeoee TEMBER 1979 CRUISE 

This cruise investigated and verified the existence of a 
satellite sensed wedge-like surface thermal feature extending 
100 km SW from Pt. Sur (Plate 1). From satellite images re- 
ceived from 18 to 26 September, it appeared that this feature 
formed off Pt. Sur on 18 September. By 20 September the apex 
of this wedge-like feature had acquired a northward, anti- 
cyclonic curvature. Between 20 and 26 September this feature 
became diffuse , a new coastal upwelling event commenced 
and a new feature was superimposed. A satellite imagery 
summary is presented in Figure 2. The best image of the 
Beateure (Plate 1) was taken at 2255 GMT, 26 September, ca-2 
hours prior to the ship's departure for the study area. 
Winds, as measured by the ship's anemometer, blew from the 
northwest averaging between 5S and 7 msec !, 

The objectives of this cruise were to obtain surface 
samples from the feature and its oceanic surroundings, to 
sample both day and night to minimize diel effects, and to 
Maxime ce time mear the cold core of the feature. A seven- 
pointed star track was planned (Fig. 3) using thermal pattern 
information derived from satellite IR imagery. Temperature, 
nitrate, and phosphate are plotted over elapsed distance in 


Figure 4. Correlation coefficients of r = 0.78 for nitrate 
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TIROS-N satellite IR image of the California coast 


for 26 September 1980. 


2 










—— 20 SEP 79 
ae 2CNSEi2 7S 
eee 02 OCT 79 


T22OW 


nn <F eae Ge ee ee ee ee ee ee ee oy + ee ei ee ee ee ee ee EPO ere re ee ee ee ee oe | a=<« eo-- ~ ew cee e ~_ SB 2 oe S Rh 


Fig. 2. Satellite features observed from 20 September to 
2 October 1979 [Johnson, 1980]. 
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Fig. 3. Track of the September 1979 cruise and outline 
(dashed line) of the oceanic front. 
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Fig. 4. Nitrate, phosphate, and sea surface temper- 
ature versus elapsed distance along the 
track of the September 1979 cruise. 
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Fig. 4 (cont'd) Nitrate, phosphate, and sea surface 
temperature versus elapsed distance along 
the track of the September 1979 cruise. 
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memonoscpniate, r= -0.91 tor nitrate to temperature, and 

r = -0.66 for phosphate to temperature were obtained (Table I). 
Nutrient ratio and temperature are plotted over elapsed dis- 
tance in Figure 5. The linear correlation coefficient between 
nutrient ratio and temperature is r = -0.50. 

The linear regression analysis of nitrate/phosphate, 
nitrate/temperature, phosphate/temperature, and nutrient 
ratio/temperature (Figs. 6, 7, 8, and 9) yielded slopes of 
momo, -2.060, -1.00, and -7.72 respectively, and x-axis in- 
tercepts of 0.45 uM phosphate, 16.68°C, 16.77°C, and 16.53°C 
mepectmvely. The x-intercept indicated that the environment 
1s nitrate limited and determined the warmest temperature for 
which these regression lines could be applied. 

The major axis of the upwelling system appeared to parallel 
and be located in the vicinity of the Sur submarine canyon 
(Figs. 10, 11, 12, and 13). The largest area of the feature 
and the part which appears to curl cyclonically is south of 
the submarine canyon. The leading edge of the curl is located 
between the Sur and Lucia submarine canyons over the shelf 
Slope, beyond the 1500 meter depth contour. 

The incremental change per kilometer of temperature and 
nutrient concentrations are plotted over elapsed distance in 
Figure 14. The oceanic front appears to be located at elapsed 
aeomamecemls slo, 39, 67, 166, 200, 312, 361, and 414 km 
(Table IT) 
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5. Nutrient ratio and sea surface temperature 
versus elapsed distance along the track of 
the September 1979 cruise. 
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> (cont'd) Nutrient ratio and sea surface 
temperature versus elapsed distance along 
the track of the September 1979 cruise. 
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Fig. 10. Sea surface temperature map for the September 1979 
cruise (contour interval, 0.5°C). Map generated 
from in situ data aided by IR imagery. 
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Surface nitrate map for the September 1979 cruise 
(contour interval, -1 uM nitrate). Map generated 
from in situ data aided by inferences from IR 
imagery. 
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Boo. 12 Surface phosphate map for the September 1979 cruise 
(contour interval, 0.2 uM phosphate). Map gener- 
ated by in situ data aided by inferences from IR 
lmagery. 
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Fig. 13 Surface nutrient ratio map for the September 1979 
cruise (contour ratio interval, 1). Map generated 
by in situ data aided by inferences from IR imagery. 
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Incremental change per kilometer of temperature, 
nitrate, and phosphate versus elapsed distance 
along the track of the September 1979 cruise. 
Frontal transit occurred at peaks. 
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Fig. 14 (cont'd) Incremental change per kilometer of 
temperature, nitrate, and phosphate versus 
elapsed distance along the track of the September 
1979 cruise. Frontal transit occurred at peaks. 
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The thermal and chemical gradients are both sharp and 
well defined despite a narrow spread of surface nutrient and 
temperature values between the ocean and the upwelling systen, 
0.01 to 4.36 uM nitrate, 0.03 to 0.96 uM phosphate, and 17.905 
to 14.80°C. The oceanic front seemed to stay in close prox- 
imity to the 15.9°C, 1.7 uM nitrate, and 0.6 uM phosphate 
isOlines. The oceanic front is outlined in Figure 3. 

Figure 15 shows the vertical thermal structure of the 
upwelling system. One hundred forty XBT probes were used 
with an average spacing of 4 km. The thermal feature is 
evident from 15 to 26, 43 to 64, 197 to 246, 372 to 383, and 
408 to 428 km. The thermocline under the upwelling system 


was continuous but upwarped to within 20m of the surface. 


B. NOVEMBER 1979 CRUISE 

The staff oceanographer at NESS began monitoring satellite 
images of the study area on 9 November. He noticed a discern- 
ible though weak thermal feature on 11 November. It extended 
35 km southwest of Pt. Sur. From satellite observation, the 
thermal feature appeared to reach its maximum seaward extent 
of at least 60 km on 20 November. By 27 November it had 
become diffuse and by 28 November was indistinguishable. Due 
to intermittent cloud cover, no satellite information pertain- 
ing to a specific feature was therefore available when the 
ship departed 28 November. A prior cruise in the area from 
26 to 27 November detected a frontal feature from in situ 


temperature measurements. 
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Fig. 15 (cont'd) Vertical temperature sections along the 
track of the September 1979 cruise [Johnson, 1980]. 
Vertical lines represent XBT drops. 
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Fig. 15 (cont'd) Vertical temperature sections along the 
track of the September 1979 cruise [Johnson, 1980]. 
Vertical lines represent XBT drops. 
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As it happened, another distinctive thermal pattern was 
seen later in images on 28 November. This feature extended 
west of Pt. Sur and grew to maximum size by 6 December. It 
advected north and apparently merged with upwelling off 
Monterey Bay around 14 December (Fig. 16). Despite inter- 
mittent cloud cover, the best satellite image was taken during 
the study period at 2310 GMT, 29 November (Plate 2). The 
average wind was 3 to 5 m-sec + from the northwest. 

The cruise plan had the ship transit south paralleling 
the coast. When the front was penetrated, an expanding square 
search would be followed to delineate the shape, center, 
orientation, and size of the feature. A ladder search pattern 
was then to be executed to provide additional transits across 
the major axis of the feature. When the expanding square 
search located a small feature south of Pt. Sur, the ship 
transited north along a front and then began a ladder search 
seaward. The first leg transected the major axis of a large 
feature located west of Pt. Sur. After the general location 
of a sharp thermal gradient was found about the feature, an 
expanding square search was begun in the vicinity of the 
leading edge (Fig. 17). 

The data herein, presented is associated with those meas- 
urements taken during this ladder search and final expanding 
square. In addition, seven oceanographic stations were 
occupied near the end of the study period. Nansen casts were 


made to sample temperature, salinity, and nitrate. 
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14 December 1979 [Johnson, 1980]. 
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Fig. 17 Track of the November 1979 cruise and outline 
(dashed line) of the oceanic front. 
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Because of contamination in the combined working reagent, 
momreliable phosphate concentrations could be measured on 
the Autoanalyzer. 

The inverse correlation between nitrate and temperature 
meae. 18 and 19) was strong, r = -0.93. The regression 
line of nitrate to temperature had a slope of -3.24 and an 
x-intercept of 14.74°C. The feature extended seaward in a 
southwesterly direction from a point near the coast ca. 20 km 
Women ot Pt. Sur. At a position ca Z0 km west of Pt. Sur, 
it appeared to curl cyclonically (Figs. 20 and 21). Unlike 
the September feature, this upwelling system did not appear 
Memexccena beyond the shelf slope to depths greater than 1200m. 

The incremental change per kilometer of temperature and 
nitrate are plotted over elapsed distance in Figure 22. The 
oceanic front appears to be transited at elapsed distances 
ieeeeetgoa, 203, 251, 257, and 270 kn. 

The thermal and nitrate gradients are sharp and well 
defined. The steepest appear to occur near the equatorward 
edge of the curl. Compared with the September cruise, the 
Menace gradient at the front is greater for the December 
feature, whereas the temperature gradients show no great 
variation (Table II). The spread of the surface temperature 
and nitrate values between the ocean and the upwelling fea- 
ture, 15.50 to 11.85°C and 0.12 to 9.92 uM nitrate, was 


greater in magnitude than the September cruise. The oceanic 
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Fig. 18 Nitrate and sea surface temperature versus elapsed 
distance along the track of the November 1979 cruise. 
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Fig. 20 Sea surface temperature map for the November 1979 
cruise (contour interval 0.5°C). Map generated 
from in situ data aided by IR imagery. 
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front seemed to closely approximate the 14.0°C and 2.5 uM 
nitrate isolines. The oceanic front is noted in Figure 17. 

Figure 23 shows the vertical profile of the upwelled 
feature based on 40 XBT profiles. The vertical thermal 
Structure is apparent from elapsed distance 179 to 216, 226 
to 247, 275 to 290, and 298 to 319 km. The thermocline under 
the upwelling system is continuous and shows no great pertur- 
bations. 

The data obtained from the seven Nansen casts appears 
in Figures 24 and 25. Because of the increase in nitrite as 
depth increases, the actual concentration of nitrate below 
the surface layer can be less than the values recorded (see 
Methods: Nutrients). 

Station 7 was taken near the major axis of the upwelling 
system, whereas station 2 sampled the warmer oceanic water 
to the south of the feature. Stations 3, 4, 5, and 6 were 
all taken in close proximity to the oceanic front. An in- 
trusion of warmer oceanic water into the central region of 
the feature appeared to have been recorded at station 1. The 
Nansen cast data showed that the upwelling system apparently 


corresponded .to a 26.4 o, tongue of water upwarping close to 


ic 
the coast and travelled seaward near the surface. The nitrate 
concentrations seem to closely agree with the density struc- 


EUre. 
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Fig. 23 (cont'd) Vertical temperature sections along the 
track of the November 1979 cruise [Johnson, 1980]. 
Vertical lines represent XBT drops; heavy lines 


represent bottom topography. 
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November 1979 cruise. 
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Ga JUNE 1980 CRUISE 

Satellite monitoring of the waters off Pt. Sur on 9 June 
indicated a pulse of upwelled water had surfaced just before 
the ship departed Monterey en route to the study area. 
Plate 3 was the best image of the surface thermal feature, 
taken at 2242 GMT 9 June 1980, four hours prior to the conm- 
mencement of the study. The wind during the study period, 
June 10 and 11, averaged between 11 and 13 m-sec + and was 
blowing from the northwest. 

Metietne Leature located by satellite thermal patterns, 
a five-pointed star track was planned so that crossings of 
the feature and coastal upwelling could be made. By halving 
the ship's speed during front transits, finer resolution data 
was to be obtained to compare the poleward and equatorward 
fronts for differences in thermal and nutrient gradients. 
The cruise plan was modified because of rough seas. After 
three crossings of the feature, the front analysis was given 
priority. A modified ladder track maneuvered the ship back 
and forth across the oceanic front. Before concluding the 
Study, the ship transected the upwelling along its major axis. 
iitessaip’s track is shown in Figure 26. 

mie=cCOhbelation Coerficients were r = 0.96 for nitrate 
to phosphate, r = -0.96 for nitrate to temperature, r = -0.92 
for phosphate to temperature, and r = -0.96 for nutrient ratio 
to temperature (Figs. 27 and 28). Linear regression analyses 


between the above parameters yielded slopes of 18.05, -4.24, 
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=0.64, and -2.07 respectively, and x-axis intercepts of 
mo7eumM phosphate, 13.47°C, 13.88°C, and 13.97°C respectively 
Wemes. 29, 30, 31, and 32). 

The major axis of the thermal and nutrient patterns are 
parallel to and directly over the axis of the Sur submarine 
canyon (Figs. 1, 33, 34, 35, and 36). The leading edge of 
the thermal feature appeared to curl cyclonically in the 
vicinity of the head of Sur submarine canyon. The crossing 
of the oceanic front appeared to have occurred at elapsed 
Gmetances 14, 36, 112, 158, 172, 206, and 222 km (Fig. 37). 

The sharp and well defined thermal and chemical gradients 
at the front seem to be less than in the November feature but 
greater than in September (Table II). However, the spread of 
temperature and nutrient concentrations is the widest for the 
June cruise; 13.77 to 9.04°C, 0.29 to 19.98 uM nitrate, and 
0.91 to 2.31 uM phosphate. The oceanic front appeared to 
closely approximate the 11.5°C, 8.8 uM nitrate, and 1.6 uM 
phosphate isolines. Figure 26 shows the oceanic front. 

The vertical temperature profile was generated from 
23 XBT profiles. The thermocline under the upwelling was 
discontinuous and appeared to surface in the vicinity of the 


oceanic front. 
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Fig. 31 Phosphate versus temperature for the June 1980 cruise. 
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Fig. 33 Sea surface temperature map for the June 1980 
cruise (contour interval, 0.5°C). Map generated 
from in situ data aided by IR imagery. 
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Fig. 34 Surface nitrate map for the June 1980 cruise 
(contour interval, 1 uM nitrate). Map 
generated by in situ data aided by inferences 
from IR imagery. 
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Fig. 36 Surface nutrient ratio map for the June 1980 cruise 
(contour ratio interval, 1). Map generated from 
in situ data aided by IR imagery. 


79 





ote 





L€G 1 LEG 2 






aetere TEMPERATURE “Clem 
cover NITRATE, pp @i/am 
scoce PHOSPHATE. »@/im 





ELraPSED O:STANCE am 


a 





L1EG 4 LEG 5 


TEMPERATURE “C/am 





| 
‘ 
NITRATE pp /am | | 
PHOSPHATE, asim | 

t 









one TEMPERATURE, “Clam 
ewoee MITAATE. 2 B/D 
eune- PHOSPHATE, 2 M/ im 





TEMPERATURE “Cam 
NITRATE, pp B/arm 
PHOSPHATE, op M/ am 





ELAPBED OISTANCE. am 


37 Incremental change per kilometer of temperature, 
nitrate and phosphate versus elapsed distance 
along the track of the June 1980 cruise. Frontal 
transits occur at peaks. 
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TABLE [I 
GRADIENTS AND OBSERVED VALUES AT THE OCEANIC FRONT 


POL LTION GH PRONT 
_ CRUIGE POULT. CRADIENT OBSERVED VALUES WHERE GrADILVts AND 
TION AT PROUT AES ROMT. VALULS MEAGURED 


‘$ 


m : 
c {4 NO - 
TC MM 3 JAM PO, 


Ee 1 
DEC 79 , 
PT suR 
1 
2 
APR 79 3 
WARM 
4. COLD 
5 
a a 
ee ay 
S 
a 
JUN 79 2 
i 
AUG 79 i 
1 3 0.6 
2 248 0.6 
3 2.3 Om 
4 1.4 0.6 
5 aS 0.3 
a 6 ee Oo 
7 0.6 0.6 
g 2.6 0.6 
9 == — 0.8 
il: 
2 
3 
NOV 79 4 
5 
6 
Me 
2 
JUN 79 3 
4 
5 
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Ve DeseUooLoN 

PeeeeoAlELLITE DETECTED THERMAL PATTERNS COMPARED 

TO IN SITU SURFACE MAPS 

Inferring the distribution of surface nutrient concen- 
trations from surface thermal patterns detected by satellites 
appear feasible from the correlation of in situ nutrients and 
temperature. In this study sea surface nutrient maps were 
generated for comparison with satellite detected thermal 
patterns formed by upwelling systems off Pt. Sur, California. 
For example, the surface nitrate and phosphate maps for the 
June 1980 feature (Figs. 34 and 35) appear similar to the 
thermal distribution (Fig. 33) with regard to the general 
location of the nutrient feature, its center, and the orien- 
tation of the major axis. However, the nutrients exhibited 
a highly structured pattern with sharp gradients within the 
feature which was not predictable from the satellite image. 
The magnitude of the difference between the distribution 
implied by the image and the actual surface distribution 
of nutrients within the feature varied among the April, Sep- 
tember, November 1979 and June 1980 upwelling systems. (The 
Other features lacked sufficient data to construct surface 
nutrient maps.) The thermal and spatial resolution of the 
AVHRR (1°C and 1.1 km respectively) could account for some 


of the difference. However, the relative ranges of the 
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nutrient values compared to the range of the temperature 
varied too much in different features at different times to 
be satisfactorily explained by the lack of resolution alone. 
B. NUTRIENT-TEMPERATURE CORRELATION AND 

"ACTIVE UPWELLING" SYSTEMS 

To better evaluate these differences, regression analyses 
were applied to in situ values of nutrients and temperature 
within the different features which formed at different times 
O@iethe year (Fig. 38). Although there were variations, the 
slopes tor all features were negative. A strong inverse 
linear correlation existed in five of the seven features 
observed since December 1978 (Table I). This implied that 
regression lines between nutrients and temperature might be 
useful in describing the distribution of nutrients with re- 
Spect to surface temperature patterns. The variations in 
Slope of these regression lines suggests that the relation- 
ship of nutrients to temperature might depend on numerous 
factors. These may include the temperature and nutrient 
Characteristics of the source water, the depth from which 
the water was advected, wind stress and bathymetric effects, 
surface dynamic processes (viz., mixing, advection, and heat 
exchange) and biological processes. Because of the complexity 
of the problem, trends were investigated with a view toward 
developing indices for prediction of nutrient distribution. 

The strong inverse correlations between temperature and 


nutrients occurring in five cruises suggested that 'active 
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38 Summary of regression lines for nitrate, phosphate, 
and nutrient ratio versus temperature 7 DEC 1978 - 
10 JUN 1980. 20 JAN, oceanic water; 7 DEC, 30 APR, 
8 AUG, and 10 JUN upwelled water from below thermo- 
cline; remaining cruises "upwelled" water from 
above thermocline. 
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upwelling' systems had occurred. This was readily apparent 
for the April 1979 and June 1980 cruises for which the ver- 
tical temperature profiles showed the thermocline surfacing 
near the edge of the surface thermal feature. This break, 
through the thermocline, would have permitted the advection 
Seenutrient-rich, "biochemically new" water [Traganza et al., 
1980], to the surface waters within the feature. This seems 
to be supported by the high nutrient concentrations, high 
nutrient ratios, and low temperatures below 10°C observed 
inside the feature (Figs. 27 and 28, respectively). 

Although satellite history of the November 1979 cruise 
Suggests the observed feature was sampled during a stage of 
"active upwelling', XBT data indicate that the system appeared 
not to have penetrated the thermocline. The thermocline under 
the feature was continuous and did not appear to be perturbed 
by the shallow (ca. 50m) upwelling system. This would have 
accounted for the moderate level of nitrate observed at the 
surface inside the feature. 

The December 1978 and June 1979 cruises did not collect 
vertical temperature data. From the wind data and the mod- 
erate levels of nitrate and phosphate, the upwelling system 
observed on the June 1979 cruise was assumed 'active' 
{[Conrad, 1980]. Because the nutrient levels were similar 
in magnitude to the November 1979 cruise, this feature could 


be representative of a shallow upwelling. 
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As the December 1978 cruise had low temperatures, high 
nutrient concentrations and nutrient ratios similar to those 
of the April 1979 and June 1980 cruises, the feature observed 
could have been associated with a discontinuity in the ther- 
mocline permitting deep upwelling. Thus, strong linear 
correlations between nutrients and temperature may be assumed 
for 'active upwelling' systems. 

C. THE RELATIONSHIP BETWEEN SOURCE WATER 

CHARACTERISTICS AND THE OCEANIC FRONT 

The ability to determine the initial condition of the 
upwelled water was considered a major factor in prediction. 
The nutrients in this subsurface source water would provide 
the maximum nutrient concentrations observed in the upwelled 
surface water. The nutrient concentrations and temperature 
were expected to be a function of depth and ambient water 
mass. The depth of the source water is most likely tied to 
the duration and strength of wind or other forcing functions, 
e.g., internal waves [Johnson, 1980]. 

Within an ‘active upwelling,’ the surface is replenished 
by vertical advection of source water. The properties char- 
acteristic of the source water are confined to the surface 
feature by the oceanic front. For every feature studied, the 
maximum gradient analysis of in situ temperature and nutrients 
(see Methods: Definition of Front and Calculation of its Grad- 
lents) showed the nutrient fronts coincided within +.6 km to 


the thermal front (Figs. 14, 22, and 37). Thus, the thermal 
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front detected in satellite imagery closely approximates 
the nutrient front. 

Since nutrients are linearly related to temperature, 
thermal gradients are well represented within a satellite 
image, and the nutrient fronts and thermal front coincide, 
the features were examined to see if the sharpness of the 
thermal front might reflect the range of nutrients within 
the feature. 

The sharpness of the thermal front was not found to be 
a dependable measure of the nutrient distribution within a 
feature. The thermal gradients at the fronts of the December 
1978, April and September 1979, and June 1980 features had 
Similar magnitude, 0.12°C/km to 0.38°C/km (Table II). Within 
the December, April, and June features high concentrations, 
eeeeto 25.1 UM nitrate and 1.0 to 2.3 uM phosphate were ob- 
served. These surface values were considered the result of 
'active' deep upwelling. However, the September feature had 
Significantly different nutrient concentrations, 0.6 to 4.3 
UM nitrate and 0.5 to 0.9 uM phosphate. The vertical temper- 
ature profile of the September cruise showed a continuous 
though greatly perturbed thermocline under the surface feature 
(Fig. 15). The low nutrient concentrations and nutrient 
ratios, together with the vertical temperature date [see 
Johnson, 1980], suggest that a shallow upwelling or mixing 
in the upper layer caused the feature. The satellite history 


of the feature revealed the thermal patterns were becoming 
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indistinguishable at the time of the cruise. The lack of 
strong correlation in the nitrate to phosphate and nutrient 
to temperature regressions also suggested that the feature 
was dissipating. The June, August, and November 1979 fea- 
tures also showed differences between nutrient distributions 


despite their similar thermal fronts. 


D. NUTRIENT DISTRIBUTION PREDICTION REQUIRES GROUND TRUTH 

Because of the inability to describe adequately and there- 
fore predict the nutrient distribution from satellite detected 
thermal patterns alone, ground truth data are needed. In 
that the linear correlation was strong between nutrients and 
temperature, an indicator of active upwelling, for the five 
features mentioned earlier, the sum effects of the initial 
conditions, boundary conditions, and dynamic processes appear 
confined to altering the slope of the nutrient to temperature 
Heoresston lime and its intercept. 
Peele FEASIBILITY OF SATELLITE PREDICTION 

AIDED BY LIMITED IN SITU SAMPLING 

A test, therefore, was conducted to determine how good a 
satellite prediction could be, using source water temperature 
and nutrient characteristics in a regression equation to 
forecast surface nutrient concentrations. The regression 
line is bounded by the coldest temperature in the feature 


and the warmest oceanic water temperature. 
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MmemuocchOumencwwccanice Water Temperature as an 
Miiecmecolmuotea + redictive Regression Line 


frewccipcragune (lj) intercepts Computed from in situ 
nutrient/temperature regression equations were compared with 
the warmest oceanic temperature to determine if a trend could 
be seen that would be useful in developing a regression equa- 
tion to forecast nutrient concentrations. The phosphate/ 
temperature T-intercept was always greater (warmer) than the 
hnitrate/temperature T-intercept. Because the oceanic en- 
vironment appeared to be nitrate-limited (nitrate goes to 
zero while phosphate is still present, see Figs. 6 and 29), 
the phosphate/temperature line would have to be extrapolated 
to reach the temperature axis, thus greater (warmer) temper- 
atures. This would be an over estimate of the actual oceanic 


water temperature. The difference, T for 


ocean” fintercept? 
active upwelling systems as defined by nutrient values and 
concentrations ranged from -0.21 to +0.61°C for nitrate/ 

temperature and -0.18 to +0.27°C for phosphate/temperature 
regression equations (Table I). Within the thermal resolu- 
tion (.5°C) of the satellite sensor, the T-intercept might 


therefore be considered equivalent to the temperature of the 


oceanic water. 


2. Hypothesis That Nutrient/Temperature 


Relation within Source Water Is Linear 
Strong linear correlations were seen between both 


nutrients and temperature and their gradients in the November 
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1979 shallow upwelling system and June 1980 deep upwelling 
ect. crom satellite history, these features appear to 
have been observed during an initial stage of development 
indicating the water upwelled from its source depth had been 
in contact with the surface for a day or less. Therefore, 
the source water nutrient and temperature characteristics 
Should be linearly related. If this is the case, and the 
source water regression line is known (from recent or his- 
torical data), this relation could be used, along with the 
aeeanhe Water temperature as a T-intercept, to provide a 
first approximation of the nutrient distribution reflected 
by the surface thermal patterns. 

In order to test this hypothesis, source water char- 
acteristics were examined using T-S and T-N diagrams from 
the shallow upwelling November 1979 Nansen casts (Fig. 25). 
The 2 and 25m T-N samples for station 7 (located near the 
coldest water in the feature) and stations 3, 4, 5, and 6 
Wiseateda mear the front) were plotted on the T-N diagram for 
Station 2, which was representative of the oceanic environment 
(Fig. 39). The 25m samples were tightly clustered on the 
Station 2 curve between the 25 and 50m depths, suggesting 
that the water had upwelled from that depth. The 2m samples 
of the ocean front stations were skewed off the station 2 
reference curve toward warmer, lower nutrient water. However, 
Station 7 varied little between the surface and 25m depth. 


Apparently station 7 was nearest the area where the source 
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Fig. 39 T-S and T-N diagram for station 2; ocean reference, 
with 25m (solid circles) and 2m (open circles) data 
from stations 3, 4, 5, 6, and 7 added to identify 
Source water. 
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water first came in contact with the surface. The ocean 
front station data were obtained in this source water after 
it had been affected by mixing, heat exchange, and other 


processes, dynamical and biological. 


fee lhe Etfects of Wind Mixing andwieat Exchange to 
Ehes Atmosphere on Surface Water near Oceanic Front 


The difference in 2 and 25m temperatures could be 
associated in part with the mixing caused by the 3 to 
sm-sec + winds during the November cruise and the atmospheric 
thermal exchange associated with a low percentage of cloud 
Bomereand Ca, 10 hours of daylight observed. Over the one 
day, the upwelled water known to have been in contact with 
the surface the longest (which was assumed to be the water 
Seeeic oceanic front) was predicted to be ca. 0.8°C warmer 
(James, 1966]. This agreed with what was observed on the 
ieemiediagram (Fig. 39). 

4. Determination of Depth From Which 

ource Water Was Advecte 

It was postulated that a nutrient source depth range 
could be determined by the depths at which the cold tempera- 
ture in the center of the feature, TO» and the temperature 
of the ocean front, Te, were found on a temperature versus 
depth profile, obtained, e.g., with an XBT. From these depths 
a range of nutrient values could be measured by in situ 
sampling. The change in nutrients versus the change in 
temperature over the depth range, corresponding to Vee 


fe ? 


a 





eauanpcea first approximation of the slope, m, for a regres- 
SHonmeequation to predict nutrient concentrations within the 
surface feature. This slope approximation could be combined 
with the oceanic water surface temperature as described 
earlier to compute the nutrient (N) intercept, b, whereby 
the total regression equation, N = mT +b, can be formulated. 
The conditions for which this approximation would be best 
Suited are those in which there were not thermal exchanges 
with the environment and changes in nutrient concentrations 
due to biological activity. Accuracy of the application to 
IR satellite data would depend on the extent of agreement 
between the surface temperature inferred from the satellite 


image and the actual surface temperature. 


weeomiinacast of Nitrate Concentration Using 
Data from the November 1979 Cruise 


The November 1979 cruise station 2 vertical tempera- 
ture profile showed the nitrate source depth range for 
mos C, Ded and 14.0°C, Te, to be 40 to 14m respectively. 
iMesnitrate range between these depths was 11.7 to 3.9 uM. 
Miecmorope tor this line, therefore, was ca. -3.7. <A second 
approximation correcting for wind mixing and thermal exchange 
with the atmosphere [James, 1966] used 13.2°C for Tp to gen- 
putea slope Of Ca. -3.12) The oceanic water surface temper- 
ature used to predict the N-intercept, b = -mT, was 15.50°C. 
From these few measurements, a predictive regression equation 


was formulated. 
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The predictive regression equation computed from 333 
in situ measurements was compared with the equations developed 
using the slope derived from source water characteristics and 
the temperature of the ocean water used as the T-intercept. 
The nitrate concentrations for 11.85° and 14.0°C were computed. 

The best fit linear regression (correlation, r= -0.93) 


based on underway in situ measurements was: 
Near 24 Tot 47.7] 4. 


where N = uM nitrate and T = °C temperature. The predictive 
regression equation, assuming no thermal exchange with the 
environment, was: 


Ny Sy 1 sata) oe SD 


The second approximation which corrected for wind mixing and 


atmospheric thermal exchange was: 


N. =P = yl Laer 0S 


a. Comparison of Predicted Nitrate Concentration 
to That Computed by Regression Analysis of 
: Observed Data 
The solution to the equations are as follows: 


for We of 11.85°C, N = 9.82 uM, N, = 13.5 uM, and N, = 11.3 uM 





it 2 
nitrate; for T, of 14.0°C, N = 2.94 uM, N, = 5.5 uM, and 
f£ 1 N. -N 
N, = 4.7 wM nitrate. The percentage of error ( x 100 


NBN 





and x 100), associated in predicting the nitrate 
concentration at the cold center or most recently upwelled 


water appears to be ca. .+37% for the first approximation, 


o 





N and ca. +15% for N, which accounted for wind mixing and 


1 Z 
thermal exchange. In the vicinity of the ocean front, the 
percentage of error increased; for Te Om a0 ae. Ny had an 
Emmon Of Ca. +98% whereas N. hadmaimchTOmrot Ca. toca. In 
both hindcasts the nutrient concentrations predicted were 
higher than those computed from the equation based on in situ 
surface measurements. Both hindcasts also deteriorated in 
reliability proceeding away from the cold center of the fea- 
ture. The predictive equation which corrected for wind mixing 
and thermal exchange with the atmosphere had significantly 
less error than the equation that did not. Because nitrate 
and temperature are nonconservative within the upper layer 

of the ocean, the 153 error derived from linear regression 
Should be considered reasonable. In all, the longer the up- 
welled water is in contact with the surface layer of the 
ocean, the more complex the interactions with the physical 

and biological environment may become, and the more unreal- 


istic the simple predictive equations. 


b. Summary of Conditions for Which This 
Test Appeared Satisfactory 


This hindcast was applied only to a shallow up- 
welling system which satellite history showed to be in an 
initial stage of development. The prediction was based on 
knowledge of the temperature of the sea surface within the 
coldest thermal pattern detected, surface temperature of the 


oceanic front detected by thermal patterns, a vertical 
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Bemperature and nutrient profile within the upper 200m of 

the ocean seaward of the oceanic front, and atmospheric 
parameters to hindcast ocean-atmosphere thermal exchange 
{[James, 1966]. In this special eee this limited 
@eound truth, it was poSsible to hindcast within 15% error 
the maximum nutrient concentration in an upwelling system 
detected by satellite imagery. It therefore seems reasonable 
Govassume that at least the major patterns of nutrient dis- 
tributions can, with further research, be inferred using 
satellite imagery and limited ground truth (i.e., data from 


buoys and AXBT's). 
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V. CONCLUSIONS 


1. Inference of nutrient distribution by satellite 
detected upwelling systems is feasible. 

2. Active upwelling systems are expected to have strong 
inverse linear correlations between nutrients and temperature. 

3. The nutrient front position can be approximated 
closely by the thermal oceanic front. 

4. The nutrient distribution within a feature can not 
be related to the sharpness of the thermal front. 

5. To predict nutrient distributions, ground truth as 
well as satellite detected thermal patterns are required. 

6. A linear regression can be used to forecast nutrient 
maxima for upwelling systems in the initial stage of develop- 
ment aided by only limited in situ data. 

7. The approximation of nutrient concentrations by linear 
regression can be improved by estimating the effects of wind 
mixing and thermal exchange with the atmosphere. 

8. With greater knowledge of source water character- 
mstres (from in Situ monitoring or historical data), stage 
of development (inferred from satellite images and in situ 
monitoring), and dynamic processes (wind mixing, advection, 
and heat transfer) a forecast of nutrient distributions with 


a surface thermal feature could be made. 
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APPENDIX A 


Listing of Surface Data: Time, Latitude, Longitude, 
Elapsed Distance, Nitrate, Phosphate, Nutrient Ratio, 
ATP, AATP/ATP, Chlorophyll, Temperature 
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